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Gravity (MOG) theory reproduces key cosmological observations without postulating the presence 



of an exotic dark matter component. MOG is a field theory based on an action principle, with 



a variable gravitational constant and a repulsive vector field with variable range. MOG yields a 



' phenomenological acceleration law that includes strong tensorial gravity partially canceled by a 

repulsive massive vector force. This acceleration law can be used to model the CMB acoustic spec- 
trum and the matter power spectrum yielding good agreement with observation. A key prediction 
of MOG is the presence of strong baryonic oscillations, which will be detectable by future surveys. 
MOG is also consistent with Type la supernova data. We also describe on-going research of the 
coupling between MOG and continuous matter, consistent with the weak equivalence principle and 
solar system observations. 

$_| 

(3JTj! Keywords: Cosmology, modified gravity, CMB acoustic spectrum, matter power spectrum, cosmic 

acceleration 

PACS: 04.20.Cv,04.50.Kd,04.80.Cc,45.20.D-,45.50.-j,98.80.-k 

> 

rt : 1. INTRODUCTION 

Why is there a need for a modified gravity theory? There is a perfectly serviceable 
model of cosmology, the so-called ACDM "concordance" model, that is not only in good 
agreement with a large body of observational evidence, it also yielded some impressive 
predictions. Nonetheless, we feel motivated to seek alternatives, in part for the following 
reasons: 

• The ACDM model requires 96% of the universe to consist of black "stuff": cold 
dark matter and dark energy, both of which may never be detectable except through 
their gravitational influence; 

• The cold dark matter paradigm runs into difficulties even closer to home, notably 
its inability to explain convincingly why the rotation curves of spiral galaxies so 
closely follow their luminosity profiles. 

The modified gravity theory we discuss here, Scalar- Tensor- Vector Gravity [1, 2] 
(STVG), has also been referred to by the acronym MOG more recently. MOG is a 
particularly interesting candidate for gravity modification in part because: 

• In the solar system or the laboratory, MOG predicts Newtonian (or Einsteinian) 
physics; 

• The MOG acceleration law is consistent with star clusters [3], galaxies [4], and 
galaxy clusters [5, 6]. 
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In the rest of this presentation, we show that MOG also appears to be consistent with 
cosmological data [7, 8]. If these results hold, MOG may prove to be a more economical 
theory (in the sense of Occam's razor) than ACDM. 

We begin with introducing MOG as a Lagrangian field theory in Section 2, also dis- 
cussing the subject of coupling between the MOG fields and matter. Next, we briefly 
introduce the phenomenology of MOG, concentrating mostly on the spherically sym- 
metric, static vacuum solution in Section 3. We then move on to cosmology: we discuss 
the MOG prediction of the acoustic spectrum of the Cosmic Microwave Background 
(CMB) in Section 4 and the galaxy-galaxy matter power spectrum in Section 5. Finally, 
we move on to the topic that is the most active area of our current research, the study of 
MOG in the presence of continuous matter such as a perfect fluid, in Section 6. We con- 
clude with a brief discussion of the most significant challenges and outlook in Section 7. 



MOG is a theory of gravity that augments Einstein's gravitational theory with a vari- 
able gravitational constant and a massive vector field with variable mass and coupling 
strength, producing a repulsive force. The theory's building blocks are: 

• The tensor field g^ v of metric gravity; 

• A scalar field G representing a variable gravitational constant; 

• A massive vector field <^ responsible for a repulsive force; 

• Another scalar field /i representing the variable mass of the vector field; 

• A further scalar field CO representing the variable coupling strength of the vector 



MOG is a theory based on a Lagrangian action principle. The MOG Lagrangian has 
three parts: the Einstein-Hilbert Lagrangian of tensor gravity, the Lagrangian of the 
massive vector field, and the Lagrangian of the three scalar fields, complete with self- 
interaction potentials: 



Here, B^ v = d^tyv — d v ^)^, and V^(<j>), Vg(G), V m (oo), and V^(ju) denote the self- 
interaction potentials associated with the vector field and the three scalar fields. The 
symbol V M is used to denote co variant differentiation with respect to the metric g^ v , 
while the symbols R, A, and g represent the Ricci-scalar, the cosmological constant, 



Although (0 is included for generality, in the solutions that we studied it turns out to be a constant. 
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and the determinant of the metric tensor, respectively. We define the Ricci tensor as 
= d«r^ v — d v r® a + r" v r^g — r"^r^ v . Our units are such that the speed of light, 

c = 1; we use the metric signature (+, — , — , — ). 

The vector field is expected to produce a repulsive force. This is not possible unless 
matter carries a vector charge. Furthermore, the vector charge must have the right 
strength to cancel out excess gravity exactly such that the effective gravitational constant 
that remains is that of Newton. This means that the coupling term must also include a 
dependence on the scalar field G. This is important for another reason as well: a scalar 
charge is required in order to ensure that the theory survives precision solar system tests 
[9]. 

We specify this coupling in the case of a massive test particle by explicitly incorpo- 
rating it into the test particle Lagrangian: 

jzf TP = - m + aG)# 5 /i w M , (2) 

where a is a function of G and q 5 is the vector charge of a test particle with mass m and 
four-velocity u^. 

This Lagrangian has been used in conjunction with the spherically symmetric, static 
vacuum solution of the MOG field equations to derive the phenomenology that we 
discuss in the next section. 



3. MOG PHENOMENOLOGY 

In MOG, the metric tensor is responsible for Einstein-like gravity, but G is generally 
greater than Newton's constant, Gm- 

The vector field is responsible for a repulsive force, canceling out part of the gravi- 
tational force; the effective gravitational constant at short range is G^. The vector field 
is massive and has limited range; beyond its range, gravity is stronger than Newton pre- 
dicts. 

The strength of G and the range ji 1 of the vector field are determined by the source 
mass. 

In the weak field, low velocity limit, the acceleration due to a spherically symmetric 
source of mass M is 

r = [1 + a - a{\ +nr)e-» r ] , (3) 

where the overdot denotes differentiation with respect to time. The values of a and \x 
are determined by the source mass M with formulas fitted using galaxy rotation and 
cosmology data: 

M f G °° A D 
a = — t= r 1 , u = — =, (4) 

(VM + E) 2 \G n J' ^ y/M 

D ~ 6250 Mq /2 kpc~ 1 , E ~ 25000 M^ 2 , = ( 1 + a) G N ~ 20G,v . 
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FIGURE 1. The acoustic spectrum of the cosmic microwave background (5-year WMAP data points 
[10] with error bars in light blue, with a moving window average also shown in purple) and the ACDM 
prediction (thick blue line). 



This acceleration law is consistent with laboratory and solar system experiments, star 
clusters, galaxies, and galaxy clusters across (at least) 15 orders of magnitude. At short 
range, fir <C 1, (3) becomes Newton's gravitational acceleration law, 

- (5) 



whereas at great distances, /ir 3> 1, we get Newtonian gravity with an "enhanced" value 
of the gravitational constant, 

r~-(l + a) — 5—. (6) 
r L 

We also used this acceleration law to investigate the MOG predictions for the cosmic 
microwave background and the galaxy-galaxy matter power spectrum, which we discuss 
below. 



4. MOG AND THE CMB 

One of the key successes of the standard model of cosmology, ACDM, is its ability 
to predict the position and size of peaks in the acoustic power spectrum of the cosmic 
microwave background (Figure 1). 

The question naturally arises: can MOG reproduce this result, especially in view of 
the fact that there exists no exotic dark matter in the MOG cosmology? 

As we were attempting to address this question, colleagues often advised us to use 
the "industry standard" cosmological code CMBFAST [11] or one of its derivatives such 
as CMBEASY [12]. When we initiated a detailed study of these programs, however, we 
found that adapting them to a variable-G cosmology is a highly nontrivial undertaking. 
A key reason is that CMBFAST uses variants of the cosmological quantity £1 (e.g., the 
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zeqpl=2 . 5d4*omegam*h*h* (2 . 7dO/tcmb) **4 
c EH (97) fitting formula for masive neutrino growth factor, 
f nu=omegan/omegam 
f cb=(omegac+omegab) /omegam 
if (fnu.gt.O.OdO) then 

apcb=0 . 2 5dO* (5 . OdO-sqrtCl. OdO+24 . OdO-f cbV) 
aktoq=(2 .7d0/tcmb)* (2 . 7d0/tcmb)/(omegam i h*h) •'h 
yf sok2=17 . 2d0'-f nu * (1 . OdO+0 . 488d0 
& *exp (-7 . OdO*l og (f nu)/6 . OdO) ) 

& *(dbl e(annunr) -aktoq/f nu)**2 

el se 

apcb=1.0d0 
yfsok2=0.0d0 
end if 

FIGURE 2. Code fragment from the CMBFAST [11] source distribution, demonstrating the difficulty 
of applying the program to a variable-G cosmology, due to the fact that D. Gp is used to represent 
constituent densities in both gravitational and non-gravitational contexts. 



baryon density to represent matter in both gravitational (e.g., structure growth) and 
nongravitational (e.g., speed-of- sound calculations) contexts. In a variable-G cosmol- 
ogy, Q. = SnGp /3H 2 may change even as the corresponding density p remains constant, 
due to changes in the value of G. Whereas gravitational relations involve the quantity 
Gp, which Q. properly represents, nongravitational relations involve p. 

This difficulty is by no means insurmountable, but it turns the adaptation of CMBFAST 
into an arduous and error-prone task. 

Other codes, such as CMBEASY, often use a version of CMBFAST as the underly- 
ing computational engine. Worse yet, the engine is often machine-translated from the 
original FORTRAN into another programming language, such as C or C++. 

It was in part for this reason that we elected to take a closer look at a promising alterna- 
tive: a semi-analytical approximation 2 developed by Mukhanov [14] that is nonetheless 
more than just a collection of fitting formulae. Mukhanov's formulation does not hide 
the underlying physics, and it becomes a relatively straightforward substitution to re- 
place, e.g., all occurrences of £1 with (G e ff/G^)n in gravitational contexts, where G e ff 
is the effective gravitational constant at the horizon. The result (Figure 3) is encouraging 
but not altogether surprising: The enhanced gravitational constant plays the same role 
as collisionless dark matter in structure growth, but dissipation is due to the baryonic 
matter density, which is the same as in ACDM. 

5. MOG AND THE MATTER POWER SPECTRUM 

Another key prediction of ACDM cosmology, confirmed by observation, is the matter 
power spectrum: the statistic of density fluctuations in the large-scale distribution of 
galaxies. 



It should be noted that the CMBFAST code base also relies on semi-analytic formulations, e.g., those 
developed by Eisenstein and Hu [13]. 
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FIGURE 3. The acoustic spectrum of the cosmic microwave background (WMAP data points with error 
bars in light blue) and the MOG prediction (thick blue line). Binned WMAP data (red) and Boomerang 
data (green) are also shown. 



This statistic can be understood using the Newtonian theory of small fluctuations in a 
self-gravitating medium. These fluctuations are governed by the equation 



d^ fC 2 s k 2 
IT 



4 + ^-4^Gp )4 = (7) 



a \ a 

for each Fourier mode 8 = 5k(/)e ik ' q (such that V 2 5 = —k 2 8; C s is the speed of sound 
in the medium). In the case of the standard model cosmology, V 2 can be obtained from 
the Poisson equation of Newtonian gravity. In MOG, the acceleration law can be used to 
derive the corresponding inhomogeneous Helmholtz equation: 

r e -tA r - f \n(r\ 

V 2 4> = AKG N p (r) + a/i 2 G N / , Zr 1 d 3 f . (8) 

J |r — r| 

The Helmholtz equation leads a shifting of the wave number in the solution to (7): 

k' 2 = k 2 + 4na 2 (^^y j - 2 . (9) 

Changes to the sound horizon scale are unaffected by the varying strength of gravity. 
However, Silk damping introduces a G 3 / 4 dependence [15]: 



*Siik = ^Silk (j^) ■ ( 10 ) 



These results can be substituted in the analytical approximations of Eisenstein and Hu 
[13], leading to the the plot shown in Figure 4. 
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FIGURE 4. The matter power spectrum. Three models are compared against five data sets: ACDM 
(dashed blue line, Q. b = 0.035, Q. c = 0.245, Q. A = 0.72, H = 71 km/s/Mpc), a baryon-only model (dotted 
green line, = 0.035, ff = 71 km/s/Mpc), and MOG (solid red line, a = 19, = 5/i Mpc -1 , £2 6 = 0.035, 
ff = 71 km/s/Mpc.) Data points are colored light blue (SDSS 2006 [16]), gold (SDSS 2004 [17]), pink 
(2dF [18]), light green (UKST [19]), and dark blue (CfA [20]). 



Comparing the result visually against the data, we can tell that the solution appears 
to have the right slope; however, unlike the ACDM prediction, MOG predicts unit 
oscillations in the power spectrum. These unit oscillations are dampened in the case 
of ACDM by the presence of matter; in the case of MOG, no such matter is present and 
oscillations are not dampened. 

Are these unit oscillations present in the data? As it turns out, it is not yet possible 
to answer this question as the resolution of the data set is not sufficient. The data are 
effectively binned using a window function. When we apply that window function to the 
MOG prediction, the unit oscillations are smoothed out, and the result shows very good 
agreement indeed with the data points (Figure 5). 

In summary, two key features of the matter power spectrum are its slope and the 
presence or absence of baryonic oscillations. MOG reproduces the right slope; however, 
unlike ACDM, MOG has unit oscillations that are not dampened by the presence of colli- 
sionless dark matter. Future galaxy surveys will unambiguously show if unit oscillations 
are present in the data. Therefore, the matter power spectrum can be key to distinguish 
modified gravity without exotic dark matter from cold dark matter theories. 

6. MOG AND CONTINUOUS MATTER 

In the preceding sections, we discussed how MOG can be used to reproduce the observed 
characteristics of the CMB acoustic spectrum and the matter power spectrum. These re- 
sults, however, were based on the MOG point particle solution (the point particle equa- 
tion of motion in the presence of a spherically symmetric, static, vacuum gravitational 
field.) 
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FIGURE 5. The effect of window functions on the power spectrum is demonstrated by applying the 
SDSS luminous red galaxy survey window functions to the MOG prediction. Baryonic oscillations are 
greatly dampened in the resulting curve (solid red line). A normalized linear ACDM estimate is also 
shown (thin blue line) for comparison. 



The question naturally arises: Is it really appropriate to use the point particle solution 
for continuous distributions of matter? To answer this question, we have to consider an 
essential feature of MOG: namely that the strength of the (Newtonian) gravitational field 
and the range of the MOG repulsive vector force both vary as functions of the source 
mass. In other words, if we combine two masses together, the resulting field will not be 
a simple sum of the individual fields, not even in a crude approximation: the theory has 
an essential nonlinearity. 

This means that we cannot blindly rely on the point particle solution to do MOG 
cosmology: it is necessary to solve the MOG field equations in the presence of contin- 
uous matter, such as a (perfect) fluid model. We have taken some tentative steps in this 
direction (see also [8]). 

How does MOG couple to continuous matter? The coupling must be consistent with 
two constraints: 

• MOG must obey the weak equivalence principle (WEP) in order to not run afoul of 
many observations, e.g., Eotvos-style experiments; 

• MOG must be compatible with precision solar system observations, expressed in 
the form of the parameters of the Parameterized Post-Newtonian (PPN) formalism, 
notably the Eddington-parameters /3 and y. 
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The two Eddington parameters /3 and 7 determine deviations from the Newtonian 
potential in post-Newtonian models: 



, 00 =l-^ + 2 ,3 (11) 
goj =0, (12) 
gjk = -(l+2^)Sjk d3) 

The Eddington-parameter /3 is identically 1 for MOG. The Eddington-parameter 7 has 
the same value as in Jordan-Brans-Dicke theory [21], which can be "cured" by introduc- 
ing a scalar charge that makes it conformally equivalent to the minimally coupled scalar 
theory [9]. 

The WEP is often interpreted as a requirement for a metric theory of gravity. This 
criteria is obviously not satisfied by MOG, as material particles are assumed to carry a 
vector charge and not move along geodesies determined by the metric. However, it is 
possible to consider a more relaxed interpretation of the WEP: a theory that satisfies the 
WEP observationally must be conformally equivalent to a metric theory of gravity. That 
is to say that there must exist a conformal transformation under which any non-minimal 
couplings between matter and non-metric gravity fields would vanish. The justification 
for this relaxed interpretation is that an observer, equipped with a classical instrument, 
would not be able to choose between conformally equivalent frames of reference: non- 
gravitational laws of physics, notably electromagnetism, are unaffected by a conformal 
transformation. 

Conformal transformations add a vector degree of freedom (the special confor- 
mal transformation, a translation preceded and followed by an inversion, x 1 ^ = (x^ 1 — 
b^x 2 ) / (2 — 2b ■ x + b 2 x 2 )) and a scalar degree of freedom (dilation, jc 7 ^ = coc^); this 
agrees with the degrees of freedom to which the matter Lagrangian is expected to cou- 
ple in the MOG theory. The metric tensor is conformally invariant up to a rescaling: 
g'Mv = a -2(j _ 2b . x + b 2 x 2 ) V v - 

These considerations about the WEP and 7 can lead to a tentative general prescrip- 
tion for the coupling between the MOG fields and matter. We anticipate that the field 
equations for a perfect fluid will contain a vector charge in the form 

(j) v u v J^ = (0 — - — T^iiy, (14) 
and a scalar charge in the form 

GJ = ~\t. (15) 

Given an equation of state, we can now write down the MOG field equations in the 
case of the FLRW metric, 

ds 2 = dt 2 -a 2 {t) [(\-kr 2 y l dr 2 + r 2 dQL 2 ] . (16) 

The equations are, after setting co = const., 
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FIGURE 6. The MOG "classical bounce". Black (solid) line is a/ao ; red (dotted) line is G/Gq ; green 
(dashed) line is jJ./jM) ; brown (dash-dot) line is (a 3 p)/(agPo). Horizontal axis is in units of 13.7 billion 
years. 
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These FLRW field equations can be solved numerically, given suitable initial condi- 
tions and some assumptions. We generally ignore the self-interaction potentials: 



= v G = y M = o. 
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FIGURE 7. The MOG "cyclical bounce" cosmology with negative Vq = const. (For legend, see Fig- 
ure 6.) 



We set the cosmological constant and the curvature constant to zero: 



We assume a simple equation of state, p = wp, and we are mainly interested in the late 
"dust" universe, w = 0. 

We use the present epoch to establish initial conditions: e.g., d/a\t=t — 
2.3 x 10~ 18 s _1 , and p|t=f — 10~ 26 kg/m 3 . The solution yields a "classical bounce", 
albeit with an age problem (Figure 6). 

To address the age problem, we can consider using a non-zero value of Vq. We find 
that a negative value produces a cyclical universe, with repeated classical bounces (Fig- 
ure 7). The conditions at the present epoch (notably, the expansion rate and approximate 
density of the universe) are repeated at later times during subsequent cycles. This allows 
for the possibility of a much older universe, one that has been through several cycles of 
shrinking and expansion since its densest state. 

What about the deceleration parameter? This parameter, defined as q = —da /a 2 , char- 
acterizes the rate at which the expansion rate slows. As it is well known, observations 
of the luminosity-distance relationship of distant Type la supernovae are inconsistent 
with a flat q = 0.5 Einstein-de Sitter universe. The data may be consistent with an empty 
universe (q = 0). It is also consistent with the ACDM universe that is dominated by dark 
energy (Q^ — 0.7) at the present epoch. 

MOG can also produce good agreement with the Type la supernova data if we assume 
the existence of a small positive value of Vq. This universe has a shallow bounce, and 
at the present epoch, its evolution is such that d/a is nearly constant. The actual age of 
the universe is, therefore, not fixed by the observed expansion rate alone; fitting to the 
supernova data yields a universe that is significantly older than the "canonical" value 
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FIGURE 8. Type la supernova luminosity-redshift data [22] and the MOG/ACDM predictions. The 
horizontal axis corresponds to the q = empty universe. The MOG result is represented by a thick (blue) 
line. Dashed (red) line is a matter-dominated Einstein de-Sitter universe with £2^ = 1, q = 0.5. Thin 
(black) line is the ACDM prediction. 



of 13.7 billion years. This model offers excellent agreement with the Type la supernova 
observations (Figure 8). 

All of these models have shortcomings, but their existence shows that MOG is capable 
of producing physically plausible models of expansion, and that a classical bounce 
occurs naturally within the theory. 

7. CONCLUSIONS AND OUTLOOK 

The research of the cosmological consequences of MOG that is presented here is on- 
going. Important cosmological results, such as the matter power spectrum and the CMB 
acoustic spectrum, have been reproduced using the MOG equations of motion in the 
gravitational field of a point source. A general prescription that describes how MOG 
couples to matter (while remaining consistent with the WEP and precision observations) 
will allow us to solve the MOG field equations in the presence of matter [23], allowing 
us to re-derive the matter power spectrum and CMB acoustic spectrum results using a 
more solid foundation. 

What can MOG tell us about the early universe? A bouncing cosmology may naturally 
avoid the horizon problem, which is a major motivation for inflationary cosmologies. It 
is unclear how MOG would address the flatness problem, other than postulating k = a 
priori. 

What about Big Bang nucleosynthesis? At short range, the MOG acceleration law is 
consistent with Newton's. However, primordial isotope ratios are sensitive to the rate of 
cosmic expansion, which is governed by long-range gravity. We do not presently know 
if MOG can be consistent with observed isotope abundances. 

Some of our results were obtained by assuming a constant value for Vq, the self- 
interaction potential for the G scalar field. This assumption is ad hoc, not justified by 
theory. Perhaps improved solutions can lead us to eliminate the need to postulate a non- 
zero Vq, or alternatively, a theoretical justification for a (near) constant self-interaction 
potential can be found. 

Notwithstanding these problems and open issues, MOG appears to be consistent with 
a range of astrophysical and cosmological phenomena. The theory yields a phenomeno- 
logical acceleration law that works well across some 15 orders of magnitude in scale. 
Furthermore, there exists an unambiguous cosmological test through which MOG can 
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be falsified: the absence of baryonic oscillations in the matter power spectrum cannot be 
explained by a MOG cosmology that lacks a cold dark matter component. 
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